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摘 要 随 着 全 球 性 的 能 源 紧缺 和 环境 问题 日 益 严 重 , 通过 充分 利用 可 再 生 能 源 和 工业 余热 资源 ， 从 而 提高 能 源 利用 效 
率 是 缓解 能 源 和 环境 问题 的 重要 方式 . 有 机 朗 肯 循环 (ORC) 是 最 有 应 用 前 景 的 低 品 位 热能 发 电 技术 之 一 . 本 文 针 对 ORC 
系统 建 并 了 结构 参数 和 系统 操作 参数 同步 优化 的 换 热 设备 多 目标 优化 模型 , 采用 R245fa 为 工 质 和 板式 换 热 器 ,以 效率 最 
大 和 比 投资 成 本 最 小 为 目标 函数 . 首先 分 析 了 单个 变量 ( 共 发 压力 、 冷凝 压力 、 过 热度 、 蒸 发 器 板 间 距 、 冷 凝 器 板 间 距 ) 对 
系统 性 能 的 影响 , 然后 选取 了 系统 的 运行 参数 ( 燕 发 压 蒸 发 压力 、 冷凝 压力 、 过 热度 ) 和 换 热 器 的 结构 参数 ( 燕 发 侨 和 冷凝 
器 的 板 长 、 板 宽 、 板 间距 ) 九 个 参数 为 决策 变量 , 利用 遗传 算法 进行 ORC 换 热 设备 结构 与 操作 参数 多 目标 同步 优化 , 获得 
多 目标 优化 的 Pareto 最 优 前 沿 及 对 应 的 最 优 系 统 运 行 参数 和 最 佳 换 热 避 结 构 参 数组 合 . 
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Heat Exchange Equipment Structure and Operation Parameter 
Multi-Objective Optimization Simultaneously of ORC 
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(Guangdong University of Technology Institute of Materials and Energy, Guangzhou 510006, China) 


Abstract The decrease in fossil energy reserves and the deterioration in environment have resulted 
in a strong interest in utilizing renewable heat sources or waste heat for power generation.The 
organic Rankine cycle (ORC) is a promising method of generating power that utilizes low-enthalpy 
renewable energy and industrial waste heat. In this paper a multi-objective model is formulated for 
the simultaneous optimization of component configurations and system operating parameters. The 
objective functions are the minimization of specific power generation cost and the maximization of 
the exergy efficiency. A waste heat driven ORC is investigated using the proposed model. R245fa is 
used as working fluid. The influences of some key parameters on the cycle performance are performed. 
Then a multi-objective optimization is performed to optimize evaporation pressure, condensation 
pressure, superheat degree, plate length, plate width and plate spacing. The exergy efficiency optimal 
scheme, economic optimal scheme and the Pare to curve are achieved using genetic algorithm. 

Key words genetic algorithm; multi-objective optimization; organic rankine; thermodynamic per- 
formance; economic appraisal 


0 引 于 


随 着 社会 的 发 展 ， 能 源 需 求 其 的 增长 ， 由 资源 
的 过 度 开发 及 不 当 使 用 造成 的 日 益 严 重 的 能 源 环 境 
问题 ， 成 为 当今 世界 人 类 需要 解决 的 重大 课题 . 通 
过 采用 热 集成 、 热 回收 等 技术 对 中 低温 余热 进行 
回收 利用 ， 可 以 提高 能 源 综合 利用 效率 ， 降 低能 源 
成 本 。 其 中 ,回收 利用 中 低温 余热 资源 对 进一步 提 
升 能 源 利 用 效率 有 巨大 的 潜力 "a。 有 机 朗 肯 循环 
(ORC) 发 电 系统 具有 结构 设备 简单 、 可 靠 、 余热 回 
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收效 率 较 高 等 特点 ,是 理想 的 中 低温 热能 回收 发 电 
技术 Ba, ORC 系统 运行 优化 和 结构 优化 可 以 提高 
ORC 的 热力 学 性 能 和 降低 成 本 . ORC 的 热力 学 性 能 
由 系统 运行 参数 和 工 质 性 质 决定 , 而 ORC 的 成 本 是 
由 组 成 系统 的 设备 结构 和 几何 参数 决定 . 因此, ORC 
系统 的 全 局 优化 应 该 是 系统 操作 参数 和 设备 的 结构 、 
几何 参数 的 综合 优化 . 

由 于 热源 温度 较 低 ， 低 温 热 源 驱动 的 ORC 系 
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统 中 超过 85% 以 上 热源 释放 的 热量 需要 换 热 器 承担 ， 
并 且 ORC 中 换 热 设备 的 损 占 总 损 的 70%~90% 左 
右 56， 且 在 ORC 的 初始 投资 成 本 中 换 热 设备 投 
资 占 整个 系统 投资 的 40%~90%[" 站 。 在 能 够 实现 换 
热 要 求 的 情况 下 ， 对 换 热 器 结构 参数 进行 优化 ,， 通 
过 常规 的 强化 传 热 方 法 ， 可 以 有 效 提高 换 热 器 的 性 
能 ， 减 小 换 热 器 所 需 面积 ， 从 而 降低 换 热 器 的 成 
本 [9.19。 然而 单纯 从 换 热 器 自身 出 发 的 传 热 过 程 强 
化 和 结构 优化 设计 ， 得 到 的 再 设计 工 况 并 不 能 保证 
系统 运行 的 综合 性 能 最 优 . 

因此 , 建立 ORC 系统 操作 参数 和 换 热 器 的 同步 
优化 模型 是 十 分 重要 的 . ORC 系统 中 的 换 热 器 运行 
工 况 受 系统 的 运行 参数 、 运 行 负 荷 、 及 冷 热源 条 件 
的 影响 . ORC 系统 与 换 热 设 备 集 成 优化 需要 考虑 在 
变 工 况 条 件 下 ORC 换 热 设备 的 优化 , 具有 优化 变量 
多 的 特点 , 采用 常规 的 线性 规划 方法 难以 进行 求解 。 
利用 遗传 算法 在 解决 非 线 性 多 目标 优化 方面 的 优势 ， 
利用 多 目标 遗传 算法 给 出 了 Pareto 非 劣 解 前 沿 6 

本 文采 用 了 板式 换 热 髓 的 ORC 系统 , 上 且 以 ORC 
系统 效率 作为 热力 学 目标 和 比 投资 成 本 作为 经 济 学 
目标 , 对 ORC 系统 的 运行 参数 (蒸发 压力 、 冷 凝 压 
力 、 过 热度 ) 及 板式 换 热 器 的 结构 参数 (蒸发 器 、 冷 
凝 器 的 板 长 、 板 宽 、 板 间距 ) 进行 同步 优化 ， 以 期 求 
解 最 佳 的 运行 参数 和 换 热 设备 的 结构 参数 . 利用 多 
目标 遗传 算法 给 出 了 Pareto 非 劣 解 前 沿 , 在 此 基础 
上 进行 决策 分 析 ， 提 供 换 热 器 的 参考 设计 工 况 ， 为 
ORC 换 热 设备 的 优化 设计 提供 参考 . 1 系统 描述 图 
1 为 纯 工 质 ORC 系统 的 Ts 图 , 整个 ORC 包括 脱 
胀 机 、 冷 凝 句 、 工 质 泵 和 燕 发 右 等 . ORC 工作 过 程 
如 图 1 中 所 示 的 1-2-3-4-5-6-7-1 过 程 , 其 中 , 1-2 为 工 
质 在 膨胀 机 中 的 膨胀 做 功 过 程 ; 2-4 为 工 质 在 冷凝 器 
中 的 定 压 放 热 过 程 ; 4-5 为 泵 对 工 质 的 加 压 过 程 ;5-7 
为 工 质 在 热 回收 节 汽 发 生 器 中 的 定 压 预 热 、 蒸 发 和 
过 热 过 程 ， 


2 ORC 系统 的 多 目标 优化 模型 


2.1 系统 热力 学 模型 
蒸发 器 中 有 机 工 质 的 理想 状态 是 定 压 蒸发 过 程 ， 
按照 温度 的 变化 情况 , 采用 分 段 方式 进行 计算 . 将 车 
发 器 分 为 预 热 段 、 营 发 展 和 过 热 段 分 段 方式 见 图 1， 
Ceva = m(hi 一 hs) (1) 
冷凝 右 中 理想 状态 下 有 机 工 质 在 冷凝 大 中 是 定 
压 冷 凝 过 程 ， 按 照 温 度 的 变化 情况 ， 采 用 分 段 方式 
进行 计算 .将 冷凝 器 分 为 预 冷 段 和 冷凝 段 ， 分 段 方 
式 见 图 1. 将 蒸发 需 分 为 预 热 段 、 蒸 发 展 和 过 热 段 分 


段 方式 见 图 1， 
人 on = m(ho 人 ha) (2) 
膨胀 机 是 ORC 系统 中 的 动力 输出 部 件 , 其 热力 
过 程 见 图 1。 
Whtur 一 71， (hi h») (3) 
工 质 泵 是 ORC 系统 中 能 量 消耗 部 件 , 其 热力 过 
程 见 图 1. 
Whump = mm: (hs — ha) (4) 
一 般 地 ， 就 参数 优化 问题 而 言 ， 优 化 模型 一 般 
包括 三 个 要 素 : 目标 函数 minf(zx)、 等 式 约束 !(z) 和 
不 等 式 约束 g(7), 如 式 (5) 所 示 : 
min fi((z), fo (2)) st 
/(Z) 一 0 (5) 
9g(Z)<0 


式 中 , 有 1(z) 和 户 (z) 分 别 为 不 同属 性 的 目标 浮 数 ， 


六 
Pre-heat Evaporation 
SS Man oe 
| | i ii\ Superheat 


的 4 Teglon 


5 : A i 2 
7 3 ! 
1 i 
“14 | 

入 aaa Cold Walet —— f : i g 
cr< | 

~ Condensation region Super cool 
region 


图 1 ORC 系统 全 一 s 
Fig. 1 T'—s diagram of ORC system 


2.1 目标 函数 

在 ORC 系统 的 优化 研究 中 , 评价 指标 是 热力 系 
统 优化 数学 模型 的 目标 函数 ， 本 文选 取 合 适 的 评价 
指标 作为 目标 函数 是 多 目标 优化 的 关键 . 

1) 热力 学 目标 

效率 :通常 的 热量 平衡 和 能 量 转换 效率 并 不 能 
反映 出 的 利用 程度 ,效率 用 来 表示 余热 资源 有 效能 
的 利用 程度 , 效率 定义 为 收益 与 支付 的 比值 : 

Wtur — Whump 


2) 经 济 性 目标 
本 文 所 研究 的 ORC 中 的 换 热 设备 均 采 用 板式 
换 热 句 . 


8 期 


比 投资 成 本 (Specific investment cost, SIC): 比 
投资 成 本 定义 为 建造 该 系统 的 总 造价 与 净 发 电量 的 
比值 ， 表 征 系统 单位 发 电量 所 需要 的 投资 成 本 ， 表 
示 为 : 

Ceom 十 Clabor 


STC = 一 一 一 7 
Whet 人 ) 


式 中 ，Ccom 为 总 设备 购置 成 本 ; Clabor 为 人 力 成 本 ， 
在 缺少 资料 的 时 ,一般 取 (0.1~0.3)Ceomi 

总 设备 购置 成 本 (Ceom): ORC 系统 的 总 设备 购 
置 成 本 主要 考虑 换 热 设备 (CHp)、 工 质 泵 (Cpump) 和 
膨胀 机 (Cutu) 等 主要 部 件 的 成 本 上 -13, 可 以 表示 为 : 


Ccom CHE 十 Cpump 斗 Ctur (8) 
CH = 10, 000 十 32440.21 (9) 
1—. 
Cpump = 422WO, a 十 1.41 ( a )| (10) 
Npup 
Ctur = 6000TT07 (11) 


其 中 ， 换 热 设 备 换 热 器 总 换 热 面积 由 蒸发 才 和 冷凝 
器 两 部 分 组 成 ,根据 热 平衡 方程 计算 面积 ， 换 热 需 
的 面积 计算 公式 如 下 ， 


A= Ceva Qeon (12) 
Ueva Teva Ucon Toon 

1 6 1 
Cava Sl 二 = 不 寺中 13 
( hao 入 | ) ( ) 

1 6 TAN 
Leon 二 | 一 二 去 寺 14 
( hevid 入 hawid ) ( ) 


2.2 约束 条 件 

该 多 目标 优化 模型 中 的 约束 条 件 包括 不 等 式 约 
束 和 等 式 约束 . 其 中 ， 各 变量 的 变化 范围 界限 、 冷 
源流 体温 度 匹 配 关系 共同 组 成 了 该 问题 的 不 等 式 约 
束 ; 模型 中 各 种 设备 的 质量 、 能 量 平衡 方程 以 及 工 
质 的 物性 方程 共同 组 成 了 该 多 目标 优化 问题 的 等 式 
约束 . 
2.3 决策 变量 

本 文 将 ORC 系统 中 板式 换 热 器 (蒸发 需 、 冷 凝 
器 ) 的 结构 参数 板 长 L、 板 宽 W.、 通道 间距 5 和 ORC 
系统 操作 参数 系统 蒸发 压力 Pva、 冷凝 压力 Peon、 循 
环 过 热度 Tsup 作为 多 目标 优化 的 决策 变量 . 其 中 , 将 
蒸发 器 和 冷凝 器 与 热源 和 冷 源 的 夹 点 做 相等 固定 处 
理 , 减少 优化 变量 , 便于 优化 求解 。 


3 案例 分 析 


对 以 压力 0.8 MPa、150°C 热 水 为 热源 ， 温 度 
15°C 冷水 为 冷 源 ，R245fa 为 工 质 的 ORC 系统 的 效 
率 和 比 投资 成 本 多 目 进行 优化 设计 . 首先， 分 别 选 


贵 奇 等 : ORC 换 热 设备 结构 与 操作 参数 多 目标 同步 优化 


取 影 响 系统 热力 学 性 能 的 系统 运行 参数 ( 燕 发 压力 
Pova、 冷凝 压力 Pon、 过 热度 Tsup) 和 影响 投资 成 本 
的 板式 换 热 融 结构 参数 (蒸发 带 板 间距 beva、 冷 族 吉 
板 间 距 bcon) 为 单独 的 决策 变量 ,进行 ORC 系统 的 
多 目标 优化 , 研究 单独 决策 变 基 对 系统 性 能 的 影响 . 
然后 , 以 系统 的 运行 参数 (va、Peon、Tsup) 和 换 热 需 
结构 参数 (蒸发 融 板 宽 em、 蒸发 器 板 长 Leva、 beva、 
冷凝 器 板 宽 帮 con、 冷凝 器 板 长 Lcon、bcon) 九 个 关键 
参数 同时 作为 决策 变量 ， 进 行 ORC 系统 的 多 目标 
优化 。 
3.1 单个 决策 变量 对 评价 目标 的 影响 
3.1.1 系统 运行 参数 对 评价 目标 的 影响 

图 2~4 分 别 给 出 了 在 其 他 参数 固定 的 条 件 下 ， 
蒸发 压力 Peva、 冷凝 压力 Pon、 过 热度 Tsup 对 系统 
性 能 系统 效率 和 系统 比 投资 成 本 的 影响 关系 。 系统 
效率 随 着 蒸发 压力 增 大 而 增 大 16%( 图 2), 随 着 冷凝 
压力 的 增 大 而 减 小 39.4%( 图 3), 受过 热度 影响 不 大 
增加 0.5%( 图 入 . 系统 比 投资 成 本 随 着 蒸发 压力 和 冷 


5800|- Poon=0.2 Tp=l Was=0.5 44 
ws=0.15 bas=0.002 Ws=0.5 a3 
5750F Zou=0.15 bs=0.002 


记 5700|- oe 
v 
局 5650 上 -40 过 
5600 上 Fa | SIC | -|39 
2 
5550| 六 一 3 
37 
14 16 18 20 22 24 26 


P,,,/MPa 
图 2 车 发 压力 对 系统 性 能 的 影响 


Fig. 2 The influence of evaporation pressure on cycle 
performance 
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图 3 冷凝 压力 对 系统 性 能 的 影响 
Fig. 3 The influence of condensation pressure on cycle 
performance 
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凝 压力 的 增 大 而 增 大 4%~54%( 图 2、3), 随 过 热度 增 
大 而 减 小 1.8%( 图 4). 
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图 4 过 热度 对 系统 性 能 的 影响 
Fig. 4 The influence of superheat degree on cycle performance 
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图 5 燕 发 器 换 热 板 板 片 间距 对 系统 性 能 的 影响 


Fig. 5 The inHuence of evaporator plate space on cycle 


3.1.2 换 热 器 结构 参数 对 评价 目标 的 影响 

由 于 换 热 器 结构 参数 对 ORC 系统 的 影响 是 通 
过 压 降 和 总 传 热 系统 这 两 个 参数 间接 作用 的 ， 因 此 
仅 以 换 热 器 的 板 片 间距 为 代表 来 分 析 换 热 设 备 结构 
参数 对 ORC 系统 性 能 的 影响 . 图 5, 6 分 别 给 出 了 
燕 发 髓 和 冷凝 右 换 热 板 板 片 间距 对 系统 性 能 的 影响 ， 
板 片 间距 对 系统 性 能 的 影响 分 析 表 明 ， 两 器 板 间 距 
增 大 系统 比 投资 成 本 增加 3%~10%, 但 对 于 循环 热 
力 性 能 改善 极 小 , 仅 有 0.26% 左 右 . 
3.2 多 目标 优化 

使 用 MATLAB(2014b) 自 带 的 基于 遗传 算法 的 


多 目标 求解 工具 箱 (GaMultiObj) 对 该 模型 进行 求 
解 . 根据 问题 的 复杂 程度 和 规模 ， 多 目标 遗传 算法 
计算 参数 确定 如 下 表 1， 下 表 2 给 出 了 各 决策 变量 
的 取 值 范围 和 多 目标 遗传 算法 进化 参数 . 
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图 6 冷凝 器 换 热 板 板 片 间距 图 对 系统 性 能 的 影响 


Fig. 6 The influence of condenser plate space on cycle 


图 7 给 出 了 ORC 的 多 目标 优化 求解 Pareto 解 . 
多 目标 求解 得 到 的 非 劣 解数 目 众 多 ,决策 者 需要 对 
系统 效率 和 系统 比 投资 成 本 的 具体 看 重 程度 . 显然 ， 
由 于 目标 函数 的 相 迟 性 ， 要 想得到 系统 效率 最 大 且 
系统 比 投资 成 本 最 小 的 解 ( 即 图 7 中 A 点 )， 是 不 
可 能 的 此 处 ,距离 A 点 距离 最 小 的 Pareto 非 劣 
解 ( 即 图 7 中 B 点 ) 通过 最 小 距离 法 求 得 . 事实 上 ， 
在 缺乏 具体 的 系统 实施 参数 的 情况 下 ，B 点 一 般 是 
决策 者 经 常 选 取 的 非 劣 解 ， 它 同时 较 好 考虑 了 系统 
效率 和 系统 比 投资 成 本 的 大 小 ， 此 时 系统 的 比 投资 
成 本 和 系统 效率 分 别 为 5588.7 €/kW 和 42.95%。 田 
外 ，C 点 为 系统 效率 最 大 43.79%, D 点 的 系统 比 投 
资 成 本 最 小 5463.6 e/kW.。 可 以 看 到 ,Pareto 前 沿 上 
的 非 劣 解 越 向 两 端 延 伸 ， 其 距离 理想 点 A 的 距离 越 
大 ， 其 实 也 就 是 越 单一 目标 的 权重 越 来 越 重 ， 不 同 
属性 指标 的 综合 评价 越 来 越 弱 . 对 应 地 , B、C、D 非 
劣 解 对 应 的 决策 变量 的 数值 解 如 表 3 所 示 . 表 3 给 
出 了 设计 参数 、 系 统 效率 最 大 、 系 统 比 投资 成 本 最 小 
和 多 目标 优化 条 件 下 系统 决策 变量 的 取 值 . 表 3 与 
图 7 相 结 合 , 很 好 地 反映 了 单 目 标 优化 和 多 目标 优 
化 的 区 别 与 改进 . 多 目标 优化 中 ， 虽 然 效 率 较 单 目 
标 优 化 是 减 小 的 , 但 是 较 系 统 比 投资 成 本 却 有 下 降 ; 


表 1 变量 的 取 值 范围 
Table 1 TablelLower and upper boundsof variables 


二 Peva/ Pda, Tsup/ Weva/ 
变量 MPa MPa °° m 
下 限 1.5 0.2 1 0.5 


上 限 2.5 0.5 5 1 


Leva/ beva/ Weon/ Leon/ bcon/ 
m m m m m 
0.15 0.002 0.5 0.15 0.002 

0.5 0.005 1 0.5 0.005 
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同 理 ， 多 目标 优化 系统 比 投资 成 本 较 单 目标 优化 是 
增 大 的 , 但 较 其 效率 却 有 提高 。 


表 2 GA 基本 参数 
Table 2 Basic parameters of GA 


参数 值 

初始 种 群 大 小 100 

最 大 进化 代数 300 

交叉 概率 0.7 
变异 概率 0.001 

最 优 前 端 个 体系 数 0.3 
适应 度 函 数 偏 差 值 10-6 


表 3 不 同 优化 目标 下 的 决策 变量 优化 结果 
Table 3 Table 3 Optimization results at different 


objective functions 


决策 变量 SIC 最 小 71 最 大 多 目标 优化 
(图 7D 点) (图 7C 点 ) (理想 点 图 7B 点 ) 
Pa 1.6338 2.4999 2.3065 
Peon 0.2002 0.2 0.2001 
Tsup 4.9485 4.8718 4.8884 
Leva 0.8251 0.629 0.766 
Weva 0.1941 0.3483 0.2615 
beva 0.0022 0.0041 0.0026 
Leceon 0.6104 0.5924 0.6 
Weon 0.3454 0.4389 0.4005 
bcon 0.0022 0.005 0.0024 
SIC/(€/kW) 5463.6 5859.9 5588.7 
Nex 38.74% 43.79% 42.95% 
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图 7 ORC 系统 多 目标 优化 Pareto 最 优 解 曲 线 
Fig. 7 Pareto curve of the ORC 


4 结 论 


ORC 运行 参数 中 燕 发 压力 和 冷凝 压力 对 系统 热 
力 性 和 经 济 性 影响 较 大 ,循环 过 热度 对 系统 性 能 的 
影响 较 小 . 板 片 间距 对 系统 性 能 的 影响 分 析 表 明 , 随 
着 板 间 距 的 减 小 , 系统 比 投资 成 本 减 小 3%~10%, 但 
对 于 使 用 R245fa 为 工 质 的 ORC 系统 而 言 , 循环 效率 
改变 极 小 仅 有 0.26%, 在 给 定 条 件 下 对 于 采用 R245fa 
为 工 质 的 ORC 系统 优化 最 大 效率 是 43.79%， 最 小 
比 投资 成 本 是 5463.6 e/kW。 通 过 多 目标 优化 得 到 
Pareto 曲线 ,用 理想 点 最 小 距离 辅助 方法 给 出 决策 


者 决策 指导 , 最 优点 B 点 效率 为 42.95%， 比 投资 成 
本 为 5588.7 €/kW. 
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